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Abstract 

Radiative penguin decays provide an indirect probe for physics beyond the Standard Model and 
contribute to the determination of the CKM matrix elements. Copious quantities of B mesons 
produced at the i?-Factories permit precision measurements of radiative penguin decays. We 
review the experimental status of the radiative penguin processes b s{d)j. 
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1 Introduction 



Radiative penguin decays are flavor-changing neutral current (FCNC) transitions that are forbidden 
in the Standard Model (SM) at the tree level, but occur at the loop level. We will focus on the 
transition with a photon in the final state which happens through penguin loops. Additional 
contributions to the electromagnetic penguin loops could arise from New Physics effects such as 
new gauge bosons, charged Higgs bosons or supersymmetric (SUSY) particles. These interfere with 
the SM processes. Depending on the sign of the interference term enhanced or depleted branching 
fractions (S's) result. Moreover, due to the presence of new weak phases; CP asymmetries that are 
small in the SM may be enhanced. This would make it possible to observe indirectly New Physics 
from the study of the radiative penguin decays. 

In addition, these decays are relevant to the determination of the CKM matrix elements. The 
measurement of the photon energy spectrum, largely insensitive to New Physics, can be used to 
improve the extraction of the CKM elements Vet and Vub from inclusive semileptonic B decays. 
Also, a direct estimate of the ratio of the CKM parameters |V((i/Vfs| can be obtained from the ratio 
of the branching fractions for the processes 6 — > ^7 and b — > ,37. 

In the following sections, the b ^ followed by the b ^ dj process will be reviewed. 



2 6 ^ 57 final states 

The SM b ^ sj branching fraction is predicted to be B{b sj) = (3.73 ± 0.3) x 10~^ P at the 
next-to-leading order (NLO). The present theoretical uncertainty of ~ 10% is dominated by the 
mass ratio of the c-quark and 6-quark and the choice of the renormalization scale. New Physics 
contributions with e.g. charged Higgs exchanges or chargino-squark loops are at the same level as 
the SM ones. In the Hamiltonian they could appear as new operators or new contributions to the 
coefficients of the SM operators. 

CP asymmetries provide another test of the SM. While small in the SM (< 1%) ^21 they can 
reach 10-50% in models beyond the SM 0. 

For the exclusive decay rate B K*^ two recent NLO calculations predict SM branching 
fractions oiB{B K*-f) = (7.ll||) x 10"^ gl and ]3{B K*j) = (7.9^1^) x lO'^ 0. The errors 
are still dominated by the uncertainties in the form factors. 

Differently from the branching fraction and CP asymmetry, the photon energy spectrum in 
6 ^ 57 is not very sensitive to New Physics processes. Moments of the photon energy spectrum 
can be used to measure the Heavy Quark Effective Theory (HQET) parameters which determine 
the 6-quark pole mass (A) and the kinetic energy (Ai) jH]. These parameters are needed to obtain 
a precision value of \Vcb\ from the b civ inclusive rate. 

Moreover, the photon spectrum can be used for the determination of V^b from B — > Xuiv. To 
lowest order in Aqcd/Mb, the B Xs^ photon energy spectrum, where Xg refers to inclusive 
strange hadronic states, is given by a convolution of the parton level 6 — s- 57 photon energy spectrum 
with the light-cone shape function of the B meson, which describes all b to light -quark transitions. 
At the same order in Aqcd/Mb, the B — > Xuiu lepton energy spectrum is given by a convolution 
of the parton level b — > uiv lepton energy spectrum with the same shape function 0. Corrections 
enter at next order in Aqcd/Mb, and these are currently the subject of active investigation |S]. 
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2.1 The Exclusive Process B K*'y 

The exclusive B K*^ modes have been studied by BABAR 9 , Belle [TO] and CLEO [Tl], where 
Belle used the highest statistics sample. Utilizing kinematic constraints resulting from a full B 
reconstruction in the B rest frame provides a substantial reduction of the qq (continuum) back- 
ground. The Belle beam-constrained ^ mass distribution is shown in Figure [2?T] for all the K* decay 
channels. 

The measured branching fractions from all the experiments and the corresponding average are 
summarized in Tabled 
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Figure 1: Belle JHl beam-constrained mass distribution for the exclusive B 
final states. 



K*j for the four K* 



Due to the large theoretical errors of 35-40%, the average branching fraction measurement is 
still consistent with the NLO SM predictions. 

The direct CP asymmetry is defined, at the quark level, as: 

r(fe s-y) - r(b s-f) 

Ac'p = = . 

r(6 ^ S7) + T{h S7) 

Table ^ summarizes the measurements of the direct CP asymmetry in S — > K*y. These are 
consistent with zero and are statistics limited. 

Isospin asymmetry, Aq-i-, is calculated by Belle using the world average value t^+It^o = 1.083 
± 0.017 [ll]- The result is: Ao+ = +0.003 ± 0.045 ± 0.018, where the first error is statistical and 

^Results for exclusive B decays are typically presented using the following kinematic variables. If {E%,p^) is the 
four-momentum of a reconstructed B candidate in the overall CM (T(45')) frame, we define 

AE* = Eg - -Bboam I 

rriES (or Mt,) = V^Wm' " Pa • 

The latter is called the energy-substituted (BABAR) or beam- constrained (Belle, CLEO) mass. Signal events peak at 
AE* near GeV and rriES (Mtc) near B meson mass; whereas continuum background lacks peaks. 
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the second systematic. It is consistent with no asymmetry, having assumed equal production of 
charged and neutral i?'s at the T(4S') mass peak. The isospin asymmetry can be used to set limits 
on Wilson coefficients jl2j . 





e(S" ^ K*"j) X 10-^ 


B{B+ K*+-i) X 10-^ 


Acp 


BABAR in] (21 fb"^) 
Belle Uni (78 fb'^) 
CLEO [Hj (9 fb~^) 


4.23 ±0.40 ±0.22 
4.09 ±0.21 ±0.19 


3.83 ±0.62 ±0.22 
4.40 ± 0.33 ± 0.24 
3.76 ±0.28 


-0.044 ±0.076 ±0.012 
-0.001 ± 0.044 ± 0.008 
±0.08 ±0.13 ±0.03 


Average 


4.18 ±0.23 


4.14 ±0.33 


-0.005 ±0.037 



Table 1: B ^ K*^ branching fraction and direct CP asymmetry measurements. 

In addition to the already established B K*j decay, there are several known resonances that 
can contribute to the Xg final state. Current measurements of higher than i^*(892) mass systems 
are from Belle ^S] and CLEO ^Ij. Note that the decay B — > (j)K'y was observed recently by 
Belle [IS] for the first time. Theoretical predictions cover a wide range; results so far are consistent 
with those from a relativistic form factor model as in Ref. |16j . 

2.2 Inclusive b ^ s'j 

Two experimental approaches have been used to measure the inclusive rate for the 6 — > 57 process. 
The "fully inclusive" method measures the high energy photon spectrum without identifying the 
hadronic system Xg. Continuum backgrounds are suppressed with event shape information, and 
then subtracted using off-resonance data. B decay backgrounds are subtracted using a generic 
Monte Carlo prediction, which is cross-checked with a 6 — > stt^ analysis. 

BABAR U7j has presented a preliminary result from a fully inclusive analysis in which the "other" 
B is leptonically tagged to almost completely suppress the continuum background. CLEO ^Hl has 
published a measurement combining several techniques to reduce the background ^. 

A "semi-inclusive" method, which measures a sum of exclusive B Xg'y decays, has been used 
by both BABAR j2Qj| and Belle |2I]. The hadronic Xs system is reconstructed by BABAR (Belle) in 
12 (16) final states with a mass range up to 2.40 (2.05) GeV. This includes about 50 % of all 6 ^ 57 
final states. Continuum and B decay backgrounds are subtracted by a fit to the beam-constrained 
B mass in the same way as in an exclusive analysis. 

Figure 121 summarizes the measurements of the b ^ sj branching fraction. The theoretical error 
is quoted as the extrapolation of the inclusive rate from the measured energy range to the full photon 
spectrum. CLEO has a lower threshold (2.0 GeV) than BABAR (2.1 GeV) and Belle (2.25 GeV). 
Presently, errors are slightly larger than the theoretical uncertainty. Computing a world average is 
complicated by the correlated systematic and theoretical errors. The fully inclusive method has a 
dominant systematic error from the B decay background subtraction. The semi-inclusive method 
has a dominant systematic error from the efficiency for reconstructing the final states, including a 
correction for missing final states that are not considered. The average branching fraction reported, 
B = (3.40 ± 0.39) X 10"^, is computed assuming the systematic errors uncor related, for simplicity. 
A first attempt to consider the correlations among the errors can be found in |22j . 

^CLEO has recently published a search for baryons in b ^ S'y events |19J . setting upper hmits on corrections to 
the B and the photon energy spectrum, which are less than half of the combined quoted statistical and systematical 
errors in Ref. 
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The present B{B ^ ^sl) measurements already provide a significant constraint on the SUSY 
parameter space. For example limits on new physics contributions to ^ Xgj have been calculated 
using the minimal supergravity model (SUGRA) [231 ^'^d charged Higgs bosons p. 

So far, only CLEO [21] has measured the direct CP asymmetry, using a technique which does 
not suppress the background coming from b ^ dj decays (which is expected to have a large CP 
asymmetry). Thus, the measured direct CP asymmetry is 0.965xAcpiB — > Xs'y)+0.02xAcpiB — > 
Xd'y) = (-0.079±0.108±0.022)x(1.0±0.03). The first error is statistical, while the second and third 
errors represent additive and multiplicative systematic uncertainties, respectively. The theoretical 
expected B{B Xdj) is used. Results are consistent with no asymmetry. 

BABAR |2flj and CLEO 18' have published a measurement of the photon energy spectrum down 
to a threshold E* > 2.1 and 2.0 GeV, respectively, where E* is measured in the B ^ and in the 
laboratory rest frame, respectively (see Figure E)). 

From the measured spectrum BABAR and CLEO have extracted the first moment in the B 
rest frame, {E^), finding (E^) = 2.35 ± 0.04 ± 0.04 GeV and (E^) = 2.346 ± 0.032 ± 0.011 GeV, 
respectively. Using expressions in the MS renormalization scheme, to order 1/M^ and order 
alPo 0, BABAR and CLEO obtain A = 0.37±0.09±0.07±0.10 GeV and A = 0.35±0.08±0.10 GeV 
from the first moment. The errors are statistical, systematic (combined in the CLEO measurement) 
and theoretical, respectively. 

Moreover, CLEO has used their measured B Xg^ photon energy spectrum to determine the 
light-cone shape function. Using this information, CLEO extracts \Vub\ = (4.08 it 0.34 it 0.44 it 
0.16 lb 0.24) X 10~^ 122], where the first two uncertainties are experimental and the last two from 
theory. 



3 b ^ d-j final states 

Both inclusive and exclusive b ^ d'y decays, which are suppressed by | Vtd/Vts P~ 1/20 with 
respect to corresponding 6 — > 57 modes, have not been seen yet. An NLO calculation, which 
includes long-distance effects of tt-quarks in the penguin loop, predicts a range of 6.0 x 10~^ < 
B{B — > Xd'y) < 2.6 x 10~^ j26] for the inclusive branching fraction. The uncertainty is dominated 
by imprecisely known CKM parameters. 

A branching fraction measurement of B ^ Xdj provides a determination of | Vtd/Vts \ with 
small theoretical uncertainties. A determination of [ Vtd/Vts \ in the exclusive modes B p{'^)l 
bears enhanced model uncertainties, since form factors are not precisely known. 

The CP asymmetry predicted in the SM for the inclusive process is foreseen between ~ 7 % 
and ~ 35 % EHl- 

Studies of the b ^ d'j decays for now focus primarily on searching for the exclusive process 
B — i- p/io^. The corresponding branching fraction is predicted to be B{B — > pj) = (I.Q^q^) x 10~^ 
[3], while the CP asymmetry is of the order of 10% |5- 

From the experimental point of view, the B p{uj)'y is more difficult than B — > K*'y because 
the backgrounds are bigger since this mode is CKM suppressed and uu, dd continuum processes are 
enhanced compared to ss continuum processes. 

■^Note that in case of a semi-inclusive analysis results can be shown in terms of the mass of the hadronic system 
Xs or, equivalently, in terms of the photon energy, E^, because = — ^^^^ in the B rest frame. 



5 



Ijlllljll 



b^sy 



ib]no- lisiak 



B +stat + syst + theo 
(3.88 + 0.36 + 0.37;J23) 10'' 

(4.3 + 0.5 + 0.8+1.3)10'' 
(3.21 + 0.43 + 0.27; J J) 10'' 
(3.36 + 0.53 + 0.42; JJJ) 10'' 

(3.11 + 0.80 + 0.72) 10'' 



(3.40 + 0.39) 10 



iliiiilii 

22.533.544.555.56 
B(10') 



'o 

CQ 150 r 
140 r 
120 r 
lOOF 
80 
60 
40 
20 




BaBar 



2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 

E., [GeV] 




Figure 2: Summary oi b ^ sj branching frac- 
tions. The shaded band shows the theoretical 
prediction described in Ref. pp. 



Figure 3: BABAR pOI (upper plot) and 
CLEO (lower plot) photon energy spectrum 
in 6 — > S7 decays. 



The smallest upper limits on the exclusive decays B p{^)7 come from BABAR which 
uses a neural network to suppress most of the continuum background. The B — > K*^ events are 
removed using particle identification to veto kaons, with a K ^ n fake rate of A multi- 

dimensional likelihood fit is made to the remaining events (see Figure^ to give 90 % C.L. upper 
limits of 1.2, 2.1 and 1.0x10^^ on p^j, p~^j and lo^, respectively. Assuming isospin symmetry, this 
gives a combined limit B{B — > p^) < 1.9 x 10^^ (90 % C.L.). Limits from Belle and CLEO can be 
found in Refs. [2H1 and respectively. 

Of particular theoretical interest is the ratio B{B pj) to B{B — s- K*'j) as most of the 
theoretical uncertainty cancels and so it can be used to determine the ratio | Vtd/Vts |. 

As described in Ref. j29j . the ratio can be written as: 

where Sp is 1/2 (1) for p^ {p^) mesons, is the ratio of the HQET form factors, AR accounts 
for possible weak annihilation and long distance contributions which appear mainly in B~^ 
p^^. Eventually these can be checked by comparing and p^'y. Using the above equation the 
constraint on the CKM elements | Vtd/Vts \ is shown in the p/rj plane in Figure [3 It is not as 
tight as the constraint from Bg/B^ mixing. However, New Physics may appear in different ways in 
penguin and mixing diagrams, so it is important to measure it in both processes. 
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dates. 

4 Conclusions and Outlook 



Figure 5: Constraint on CKM unitarity triangle 
obtained in Ref. [21] from the upper limit on 



A review of recent experimental results of radiative penguin decays b s{d)'j is presented. The 
6 — > S7 process in the exclusive and exclusive final states is well established and more statistics 
can be used to improve the limits (or indirectly find evidence) of New Physics and to improve the 
measurement of the photon energy spectrum, toward a better determination of the CKM matrix 
elements. 

There is not yet evidence of 6 — > ^7 decays but BABAR and Belle expect to collect ^ 500 fb~^ 
by 2005. This should be sufficient to observe B — > p'j. It may also be feasible to measure the 
inclusive b ^ dj rate. For the measurement of |Vtd/Vfs|, the ratio of b ^ dj to b ^ has much 
smaller theoretical uncertainties than the ratio of the exclusive decays. 
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